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G\ . Abstract 

in 

I/"") \ The LHC has started to constrain super symmetry-breaking parameters by setting 

| bounds on possible colored particles at the weak scale. Moreover, constraints from Higgs 

physics, flavor physics, the anomalous magnetic moment of the muon, as well as from 
searches at LEP and the Tevatron have set additional bounds on these parameters. Renor- 
^ , malization Group Invariants (RGIs) provide a very useful way of representing the allowed 

parameter space by making direct connection with the values of these parameters at the 
messenger scale. Using a general approach, based on the pMSSM parametrization of the 
soft supersymmetry-breaking parameters, we analyze the current experimental constraints 
to determine the probability distributions for the RGIs. As examples of their application, 
we use these distributions to analyze the question of Gaugino Mass Unification and to 
probabilistically determine the parameters of General and Minimal Gauge Mediation with 
arbitrary Higgs mass parameters at the Messenger Scale. 



*http: / /theory.fnal.gov 



1 Introduction 



The Standard Model (SM) provides an excellent description of all experimentally measured 
observables at present. Mass generation relies on the Higgs mechanism, which is based on 
the introduction of an elementary scalar field transforming in the fundamental representation 
of the SU{2)i J group. The vacuum expectation value (vev) of this scalar field sets the weak 
scale, which is then proportional to the magnitude of the square root of the negative squared 



mass parameter in the scalar Higgs potential [Jj 



2|j. The SM provides no explanation for the 



magnitude of this mass parameter, which is sensitive via radiative corrections to new physics 
at high scales. 

The Minimal Supersymmetric Extension of the Standard Model (MSSM) has most of the 
virtues of the SM (sj-jsj. Apart from a loop factor, the magnitude of the Higgs mass parameter 
is determined by the size of the supersymmetry-breaking parameters of the third generation 
squarks. These also determine the value of the SM-like Higgs mass at the loop level. Values of 
the third generation squark masses of about 1 TeV lead to SM-like Higgs masses in the 115— 
130 GeV range jf]]- [l5|. Hence, recent hints of a Higgs mass of about 125 GeV are consistent 
with MSSM predictions |l6 |. 

The supersymmetry-breaking mass parameters depend on the unknown mechanism of super- 
symmetry breaking and on the messenger scale, at which supersymmetry breaking is transmit- 
ted to the observable sector. Recent experimental bounds from the LHC set strong constraints 
on colored particles at the TeV scale, and therefore on the parameters of minimal models of 
supersymmetry breaking. 

Several works have studied the relationship of the supersymmetric mass parameters between 
the messenger scale and the weak scale 12J- 32J. It would be very useful to have a method that 
allowed us to set bounds on the supersymmetry-breaking parameters at the messenger scale, 
independent of the unknown supersymmetry-breaking scheme and of the unknown value of the 
messenger scale. Renormalization Group Invariants (RGIs) 33[- 41[ provide such a method. 
Determination of the value of the RGIs at the TeV scale sets their values at the messenger scale. 
One can then use the information provided by the RGIs to set constraints on general classes 
of models 39[, 4l|. An exhaustive analysis of the RGIs for different supersymmetry-breaking 



scenarios is performed in Ref. 42]. 

The effects of various preLHC and LHC results on thephenomenological MSSM (pMSSM) 
parameter space have been studied in detail in Refs. 2jJ-[27|. In this article, we use the pMSSM 
parametrization of the soft supersymmetry-breaking parameters 2l| to determine the current 
probability distribution of the RGIs at the TeV scale. We shall compare the situation before 
and after constraints from the LHC are imposed. 

To illustrate the power of this framework, we will use the pMSSM RGI probability distri- 
butions to analyze three particular issues: 
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• Possible scale of Gaugino Mass Unification. 



• Messenger scale parameters in a realization of General Gauge Mediation. 

• Messenger scale parameters associated with Minimal Gauge Mediation. 

The probabilistic interpretation of the RGIs can be applied to other quantities of interest in 



the MSSM using for example the analysis presented in Ref. [42] . 

In section 2 we list the RGIs to be used in this paper, outlining the methodology to be used in 
our analyses. We then compute the RGI probability distributions obtained by imposing current 
experimental constraints. In section 3 we study the question of Gaugino Mass Unification and 
the consistency of the scale of this Gaugino Mass Unification with experimental constraints. 
In section 4 we look at General Gauge Mediation and determine the probability distribution of 
the relevant parameters of this model. Section 5 discusses the probability distributions for the 
Minimal Gauge Mediated parameters. We reserve section 6 for our conclusions. Details about 
our probability analysis are given in Appendix A. Appendix B gives the specific definition of the 
pMSSM. Appendix C gives the inverted relationships between the soft masses of the pMSSM 
and the RGIs. Appendix D lists these in the case of flavor-blind models. 

2 RG Invariants: Probabilistic Interpretation 

2.1 RGI-pMSSM Basis 



There are 14 relevant RGIs, analyzed in Ref. 39(,4l|], involving the soft supersymmetry 



breaking parameters, which we will use as the basis of our current work. These are summarized 



at one-loop accuracy in Table [U two-loop corrections were studied in Ref. 39| and shown to 
be of order of a few percent or less. Moreover, there are 2 RGIs relating only the gauge 
couplings (I g2 and I g3 ), which we can use to redefine the other 12 RGIs in terms of just the 
soft masses and the scale. These soft masses, ignoring possible small flavor dependence of the 
sfermion and Higgs mass parameters, are given by a total of 17 scalar masses plus 3 gaugino 
masses. One can make the additional well-motivated assumption of degeneracy for the first 
and second generation sfermion mass parameters. In such a case, one is left with 12 scalar 
masses. Therefore, the 12 RGIs, which are linearly independent, can be inverted to give 12 soft 
supersymmetry-breaking masses in terms of these RGIs as a function of 3 given soft masses. 

In the pMSSM, apart from the 15 soft supersymmetry-breaking parameters discussed above, 
there are 3 soft supersymmetry-breaking parameters, Af (/ = t,b,r), denoting the mixing of 
the left- and right-handed third generation sfermions, and tan /3, the ratio of the Higgs vevs (or 
equivalently the soft supersymmetry-breaking parameter B^). These 19 parameters then give 
the complete basis for the pMSSM. 



2 



Table I: 1-Loop RG Invariants in the MSSM 
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The Higgs soft supersymmetry-breaking parameters may be determined, up to loop correc- 
tions, as a function of the Higgsino mass parameter, /x, the CP-odd Higgs mass, ttia, and tan (3. 
The tree-level expressions for \i and m^ in terms of m# u and m# d are given by: 

„ mi, tan 2 8 — 171% 1 „ 
" = V»*fl 2"*' <2 ' 1) 

2 (m^ - m 2 ^) 
ml. = 2 ■ (2.2) 

A cos2^ z V ; 

These can be inverted to give the Higgs soft supersymmetry-breaking mass parameters: 

1 



m 2 



2 Hu = - [(1 + cos 2/3)m 2 4 + m| cos 28 - 2/i 2 ] , (2.3) 



m 2 Hd = - [(1 -cos2/3)m 2 1 -m|cos2/3-2/i 2 ] . (2.4) 

Using the above relations, one can define a 1-to-l correspondence between the 19 pMSSM 
parameters and the 12 RGIs, complemented by /i, m^, tan/3, the 3 mixing parameters A t , Ab 
and A T and one third generation squark mass parameter, for instance mQ 3 . The expressions 
for the soft masses in terms of the RGIs are given in Appendix C. 



2.2 Methodology 



In Ref. [25j, the probability distributions of the 19 pMSSM parameters were computed, an- 
alyzing the differences between the results for these distributions considering some preLHC 
measurements (listed in Table [TT]) and after including various 1 fb _1 LHC results (listed in 
Table [TTT]) . We refer the reader to Ref. 25j for specific details about the likelihood analysis. We 



shall use the set of pMSSM points and their corresponding preLHC and LHC likelihoods from 
Ref. 25j and obtain the probability distributions for the RGIs in Table [I] projected to 5 fb _1 of 



LHC data. 

In Ref. [25[, a flat prior for the all the soft supersymmetry-breaking parameters was used. 
The scalar masses were varied between and 3 TeV. The gaugino masses and the \i parameter 
were scanned between -3 and 3 TeV, and the mixing parameters, Af, were scanned from -7 to 
7 TeV. The range of tan /3 considered was 2 to 60. 

The RGIs, are functions of the soft mass parameters and therefore a flat prior for the later 
does not imply a flat prior for the RGIs. In particular, even in the simplest cases, the fact that 
the parameters have boundary values imply that certain regions are preferred. As an example, 
consider the subtraction of two mass-squared parameters, 

f(a,b)=m 2 a -m 2 b . (2.5) 
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Table II: The pre-LHC experimental results that are the basis of our pMSSM parameter scan using 
Markov Chain Monte Carlo (MCMC) sampling. We re-weight a posteriori with the limit 
BR(B S ->■ /i/i) < 1.08 x 10~ 8 at 95% CL 0. However, this hardly has any effect. In 
evaluating the Higgs mass limit, we apply a Gauss-distributed theoretical uncertainty with 
a = 1.5 GeV to the mn computed by SoftSUSY. 
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BR(b -»■ s 7 ) 


(3.55 ±0.34) x 10" 4 [44, 45, 


Gaussian 


o 
Z 
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1/(1 + exp( 001D J) 


3 




1.66 ± 0.54 [46] 


Gaussian 


4 


A „ 


(28.7 ± 8.0) x 10 [e^e J [47] 
(19.5 ±8.3) x 10- 10 [r+r-] [47] 


Weighted Gaussian average 


5 


m t 


173.3 ± 1.1 GeV [48] 


Gaussian 


6 


m b (m b ) 


4.19^;- GeV[46] 


Two-sided Gaussian 


7 


a,(M z ) 


0.1176 ±0.002 [49] 


Gaussian 






LEP&Tevatron 


rrih sampled from Gauss(rrih, 1.5) 


8 


m/, 


(HiggsBounds [50]) 


Lg = 1 if allowed. 
Lg = 10~ 9 if excluded. 


9 


sparticle 


LEP 


L 9 = 1 if allowed 




masses 


Neutral LSP (MicrOMEGAs [51]) 


L 9 = 10~ 9 if excluded 



If both m a and m& have a uniform distribution between and 3 TeV, it is clear that the 
probability of f(a,b) ~ ±(3 TeV) 2 will be much smaller than the probability of f(a,b) ~ 0. 
This is because in the former case one of the mass parameters has to be equal to while the 
other is 3 TeV, while the later situation comprises of all cases in which m a ~ m&, independent of 
their value. Therefore, in order to determine the probability distributions of the invariants (and 
other functions considered later), one should normalize them such that they can be compared 
to a flat prior for the functions under consideration and not for the masses. In order to do 
this, we have re-scaled the experimentally weighted probability distributions of the RGIs by 
the probability distributions for these quantities obtained by varying the mass parameters 
with a flat uniform distribution in the region originally scanned. We will refer to the later 
distributions as the "Flat Un-weighted" distributions. The details of the exact procedure are 
given in Appendix A. 

The results are shown in Figs. [1] and [2j The shaded green region represents the flat un- 
weighted distribution for the RGI being considered. For the Lb,s, which depend only on, and 
are linearly proportional to the gaugino masses, this distribution is flat (apart from a small 
variation with the gauge couplings). However, for the other RGIs, these distributions acquire 
definite features. The green line represents the probability distribution obtained considering 
only the pre-LHC constraints listed in Table [XXJ These depend heavily on the constraints on 
the weak eigenstates coming from LEP and (g^ — 2), the bounds on the CP-odd and charged 
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Table III: LHC measurements used in the current study. The ax variable is effective in suppressing 
background from light-quark QCD. SS 21, and OS It denote same-sign and opposite-sign 
dileptons, respectively. The olt \E^\, SS [53|], and OS 54] results were published by the 
CMS Collaboration. We re-weight a posteriori with the limit BR(B S — > fi/j,) < 4.5 x 10 -9 



at 95% CL [55(]. This has an effect only on Dz, which depends on the Higgs mass parameter 



m lf d - We also update the Higgs bounds, imposing the constraints from the CMS di-photon 
searches 561 ] . which do not have a strong effect on the probability distribution of the soft 
supersymmetry-breaking parameters. 



j 


Analysis and search region 
(values in GeV) 


Observed event count 
(N 3 ) 


Data-driven SM 
BG estimate 
(Bj ± SBj) 


1 

2 
3 
4 
5 
6 
7 
8 


ttT hadronic, 275 < Ht < 325 
ttT hadronic, 325 < Ht < 375 
ax hadronic, 375 < Ht < 475 
ax hadronic, 475 < Ht < 575 
ay hadronic, 575 < Ht < 675 
ot hadronic, 675 < Ht < 775 
ot hadronic, 775 < Ht < 875 
a-r hadronic, 875 < Ht 


782 
321 
196 
62 

21 

6 

3 

1 


787.4t^;g 
310.4^ 2 4 4 
202.1^-e 

60.4±|;g 
20.3ii? 

7 7+0.8 
' • ' -0.5 
q O+0.4 
O.Z_ Q 2 

9 O+0.4 


9 


SS 2£, H T > 400, $ T > 12 


1 


2.3 ± 1.2 


10 


OS 2£, H T > 300, $ T > 275 


8 


4.2 ± 1.3 




Observable 


Experimental result 


Likelihood function 


11 


BR(B S /i/i) 


< 4.5 x 10" 9 [55, 


V(l + exp(f^)) 


12 


m h 


m H S*V } [56] 


L12 = 1 if allowed. 
L 12 = 10~ 9 if excluded. 



Higgs masses and third generation masses coming from the BR(b — > S7) and the LEP/Tevatron 
Higgs results. The red line, instead, represents the probability distributions obtained after the 
LHC results are considered (Table HTT1) . 

The details on the computation of the final resultant distribution we label as "p^lExp.) 
Reweighted" are given explicitly in Appendix A. However, the method can be simply understood 
by noting that, as discussed in Appendix A, the ratio of the difference of any 2 probabilities 
with the flat distribution, (pi — pf)/(p2 — p ), is preserved when the scan range on the original 
masses is changed. Therefore, we first subtract the probability distribution determined after 
the LHC measurements (red line) from the one obtained with the flat masses prior (shaded 
green) and then shift this distribution by the minimum, setting the minimum at zero. Once 
this is done, we re-scale the distribution such that the total probability is 1. This then gives 
the dashed black line denoting our final resultant distribution. This distribution is flat in 
regions not scanned or impacted by experiment and enhances and reflects the actual effect of 
the experiments on the RGIs. Larger values of this distribution highlight the regions where 
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experimental input has increased the likelihood and values less than the flat probability show 
regions where experiments disfavor model space. 

Due to this method of re-scaling and subtracting the probabilities, one has to be careful 
when using these distributions to calculate resultant quantities of interest (as will be done, for 
example, when calculating product probabilities). One needs to convert the distribution back 
to a true probability via the scale factor SF (labeled "SF" in plots) defined in Appendix A: 

p(#|post-LHC) = p(8\Flat) + — \p(6\Exp.) Reweighted - p(6\Flat)] , (2.6) 

SF 

where p(9\Flat) = l/(no. bins). Unless otherwise noted, no. bins = 100 in all plots, implying 
p(6\Flat) = 0.01. 

2.3 Results 

Here we will discuss the probability distributions presented in Figs. [T] and |5] for the different 
RG Invariants. We don't analyze Dy 13H and Dy a since these two depend on almost all the soft 
masses and the current experimental bounds on these combinations are too weak to show an 
effect on the probability distributions. 

The first three distributions displayed in Fig. Q] are for the three 1^ which are equal to 
the gaugino masses divided by the square of the corresponding gauge coupling, Mi/gf. The 
LEP constraints on charginos, sleptons and gluinos, together with the requirement of a neutral 
particle to be the lightest supersymmetric one, lead to a preference towards low values of the 
bino mass, Mi, increasing the probability for small Ib 1 - The LHC modifies this distribution 
indirectly, through the updated bounds on the gluino and squark masses. This is due to the 
requirement of having a neutral particle as the LSP: When at least one of the squarks and/or 
the gluino is light, one neutralino or a sneutrino is forced to be even lighter. Of all the neutral 
particles, the bino is the only one that is not related to the mass of other charged particles 
and therefore is not pushed to larger values. Hence, the bino can be very light increasing 
its probability of being the LSP. For heavier squarks and gluinos, the neutral particles can 
be heavier and consequently the bino mass probability distribution moves to larger values. 
Regarding Ib 2 , the LEP constraints on chargino masses, together with the bounds on 
restrict small values of M 2 , while leading to a preference for values of M 2 near the weak 
scale. The LHC does not significantly modify this constraint. Since we have not implemented 
the Tevatron bounds on the squark and gluino masses, the pre-LHC constraints on Ig 3 are 
dominated by the indirect effect of requiring that a gluino cannot be the LSP, which therefore 
disfavor the region of small values of this quantity. The LHC SUSY searches further constrain 
values of M3 up to ~1 TeV, which is clearly shown in the Ib 3 distribution. 

Looking at 1m x , also displayed in Fig. [TJ, the lower bound on the slepton and squark masses 
from LEP disfavor the lowest values of this RGI. The LHC further strengthens this trend by 
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Figure 1: Distribution of the RGIs before and after the LHC constraints are added (green and red 
lines), flat distribution (shaded green) and subtracted probability distribution (dashed 
black line). 
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Figure 2: Distribution of the RGIs before and after the LHC constraints are added (green and red 
lines), flat distribution (shaded green) and subtracted probability distribution (dashed 
black line). 
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increasing the bounds on the squark masses. Im 2 is strongly shifted to lower values by 2), 
which leads to a preference towards small values of M2 and the left-handed slepton masses. The 
LHC data does not modify this trend in any significant way. Finally, for Im 3 , the previously 
discussed bounds on the gluino (and similar bounds on the squark) masses, push them to larger 
values, leading to a preference towards non-zero values of this RGI. The LHC, with significantly 
larger bounds on the gluino mass, leads to a further preference towards larger values of this 
RGI. The asymmetry between positive and negative values comes from the fact that the gluino 
has a larger cross section and is constrained to be heavy even in the case in which only the 
third generation squarks are lighter than the gluinos. Individual squark species, on the other 
hand have lower cross sections and therefore have a higher probability of being lighter in the 
pMSSM 0. 

In order to understand the behavior of the other four invariants displayed in Fig. El it is 
convenient to analyze the results of Ref . j^J . The bound on the Higgs mass leads to a preference 
for larger values of the third generation masses, beyond the bounds on the first and second 
generation masses obtained at the Tevatron. Since the Higgs mass bounds are approximately 
symmetric in their dependence on mg 3 and m U3 and the negative weight of rriQ 3 on Db 13 is twice 
as large as the positive one of m U3 , we see a preference towards negative values of D Bl3 . The 
LHC SUSY searches have not yet changed this tendency in a significant way. Regarding Dl 13 , 
the preference towards small values of (g^—2) lead to a preference towards small left-handed 
second generation sleptons, beyond the LEP constraints, and therefore towards lower values of 
this RGI. The LHC, again, does not have a strong impact on this distribution. Lower values of 
the left-handed slepton masses also affect the pre-LHC distribution of D Xl , pushing it to lower 
values. At the LHC, there is a somewhat stronger constraint on the left-handed squarks with 
respect to the right-handed ones, leading to slightly lower values of D Xl . Finally, Dz is strongly 
dominated by the bounds on the CP-odd Higgs mass coming from B s — > fifi, which push rn 2 Hd 
to larger values and Dz to lower ones. 

3 Gaugino Mass Unification 

Gaugino Mass unification is a common feature of models in which supersymmetry breaking 
occurs at scales larger than the GUT scale. In such a case, up to threshold corrections, one 
should expect that due to the extended gauge structure, the gaugino masses associated with 
the SU{3) C , SU{2) L and U{1) Y unify at the GUT scale. 

At scales lower than the GUT scale, however, threshold corrections can be large and could 
lead to quite different values of the three gaugino masses. This happens, in particular, if the 
gaugino masses receive large contributions induced by gravitational interactions governed by 
the scale anomaly. These contributions are proportional to the /3 function coefficients of the 
respective gauge couplings. 
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Further contributions to the gaugino masses may come from gauge mediation, induced by 
messengers charged under the SM gauge groups and with direct coupling to the supersymmetry- 
breaking sector. In minimal models, the gauge-mediated gaugino mass contributions at the 
messenger mass scale are proportional to the value of the square of the gauge couplings at the 
same scale. 

The simplest contributions to the gaugino masses at the messenger scale are then given by 
Mi = A ft g] + B g 2 8(M 2 mcss - Q 2 ) + ^ g 2 (3.7) 

fi'GUT 

where the coefficients A, B and Mi/ 2 parameterize the anomaly-mediated, minimal gauge- 
mediated and minimal SUGRA-mediated contributions to the gaugino masses. We have in- 
serted a 9 function to denote the fact that the gauge mediated contribution is only relevant at 
energy scales, Q, below the messenger mass scale, M mess . 

The condition of gaugino mass unification can be written in terms of RGIs. Indeed, assuming 
that the gaugino masses unify at some scale, M unif , 

02 01 ~ M 9 ' 

~7 — ~7 
91 92 

P3 'T P ;! B ' = M„ , (3.8) 

91 9i 

where 0i = {33/5,1,-3} for i = {1,2,3}, gi are the gauge couplings at the gaugino mass 
unification scale and M g is the common gaugino mass. The denominators in the above equation 
are nothing but I g2 and — 3I g3 , respectively. The value of the gauge couplings at the gaugino 
unification scale may be obtained by just dividing the above expressions by the corresponding 



invariant. In particular 39] 



Si 2 (M, 



unif J 



Pi — PiIb 2 /Ibi ^ 02 — PiIb 2 / Ib-l 
T g% 2(/3 2 -A) 



9l(M nmi ) = - ^;/ R J i\ Bl (3-9) 

— 61 g- 3 Z \P3 — Pi) 

where we have used the fact that g%uT — 1/2. From the equality of the first and second line in 
Eq. f)3.9p . one can see that gaugino mass unification requires that 

(A» - 02)I Bl - {03 ~ 0i)I B 2 + (02 ~ 0i)Ib 3 = 0, (3.10) 



or, inserting the numerical values of the 0{ coefficients [39(, 



12I B2 - 5I Bl - 7I B3 = 0. (3.11] 
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Using the expression for the gaugino mass, Eq. 13. 7\ we get that at the weak scale, 

I Bi = Af3i+ (B + M 1/2 /g 2 GXJT ) =A(3 l + C (3.12) 

where we have joined the scale and gauge mediated contributions, C = (B + Mi/ 2 /# GUT ), smce 
they cannot be distinguished at low energies. Note that this is similar to the case of Mirage 



57| [58[ [59|. Observe that the condition given in Eq. (I3.10p is always satisfied in 



mediation 
this case. 

Interestingly enough, inserting Eq. 13.121 in both the expressions for gf(M un i{) in Eq. (13. 9p . 
the same equation is obtained, giving the necessary condition for the unification of gaugino 
masses at some scale, 

In addition, the above expression is independent of (3 2 ,3- 

In general, for positive values of A and C, depending on which supersymmetry-breaking 
mechanism is dominant, the apparent gaugino unification scale can vary from the infrared to 
the GUT scale. In order for gaugino mass unification to take place at a physical scale, however, 
we need that 0.5 > g\{M un \{) > 0.2, which sets interesting constraints on the values of A and 
C. For A = 0, one gets that unification is at the GUT scale. For C — 0, instead, one obtains 
that unification occurs for gf(M uni {) ~ 0, which is an un-physical value, and for which Ig 1 
diverges unless the gaugino mass also vanishes. Let us remark, however, that the unification 
of gaugino masses obtained by extrapolating the RG evolution into un-physical scale values, 
at which the physical spectrum is not the MSSM one, could still say something relevant about 
the supersymmetry-breaking mechanism. For example, the unphysical gf(M un a) ~ for C = 
is the expected apparent unification value in anomaly mediation scenarios. It is therefore very 
interesting to use the above expressions, Eq. (I3.9p . to check the consistency of gaugino mass 
unification, and to obtain information about the scale at which it may occur. 

Let us elaborate further on the above point. Although we computed the gaugino unification 
scale by using gf, we could have used any other gauge coupling, and the condition in Eq. (13.131) 
would be the same, with g\ and /3i replaced by the corresponding gf and The fact that the 
conditions one obtains are consistent with each other can be obtained by rewriting Eq. (13.131) 
for any gf, in the following way 

2fl+4ftV ( 3 - 14 ) 



<?f(M miif ) V C 

This has the correct form for the evolution equation for the gauge couplings from the scale 



12 




Figure 3: Distributions of the value of the gauge coupling g\ at the Gaugino Mass unification scale 
before and after LHC constraints (green and red lines), flat distribution (shaded green) 
and subtracted probability distribution (dashed black line), no. bins = 50 for these plots, 
implying p(9\Flat) = 0.02. 

Mqut, where 1/gf = 2, to other energies, provided we interpret 

M GUT \ 2 A 



or, equivalently 



M unif = Mqut exp I -16tt 2 ' 



C 

M GUT exp 



57T 



(3.16) 



Since (3^ is negative, for large values of A/C, the effective scale, Q un a, may be below the QCD 
Landau pole and therefore <j$, from Eq. (13.14p becomes negative, and so clearly un-physical. 
Negative values of g\ 2 may also be obtained for negative values of A or C. 

In order to analyze the experimental impact on the scale of Gaugino Mass Unification, we 
have studied the two possible independent determinations of g\{M un \i) coming from the ratios 
Ib 2 /Ib! and Ib 3 /Ib!, Eq- (13.91) . respectively. If gaugino masses unify at a certain scale, those 
two determinations should lead to the same value of gf{M nn ^). Since the unification scale is 
not necessarily the messenger scale, we have only restricted the value of the gauge coupling, 
g^(M unif ), to lie between and 1. 

Fig. [3] represents the probability distributions of gf(M nni f) obtained by the two ways de- 
scribed above (Eq. 13. 9p . The green, red and black curves and the shaded green area have the 
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p(0 | post-LHC), SF = 1 

M unif (GeV) 




9?( M j-( i/| g2 )(i - (33/5)(| B /y) 



p(G | pre-post LHC), SF = 1 .03 




Figure 4: Probability Distributions for the gauge coupling, gf, at the Gaugino Mass unification scale. 

Left: After LHC constraints. Right: Difference between pre and post LHC probabilities. 
There are 50 bins for each axis in these plots, implying p(9\Flat) = 0.4 x 10 -3 . 

same interpretation as the one in Figs. [T] and [2J The results are very interesting, since values of 
the gauge coupling of about its weak scale value gf(M uni f) ~ 0.2, are clearly disfavored, while 
unification at the GUT scale gf(M uni {) ~ 0.5 or at values consistent with anomaly mediation 
(?i(M uni f) ~ are somewhat preferred. 

Fig. H] shows a two-dimensional representation of these results, comparing the results for 
g\ (M un if) obtained by the two different equations. The left panel shows the results after LHC 
constraints are used and the right panel shows the difference between pre- and post-LHC. 
Dark (light) blue is strongly (weakly) disfavored, while the green regions provide an acceptable 
description. Yellow (red ) is weakly (strongly) preferred. The black diagonal line shows the 
(equal) values that should be obtained for Gaugino Mass Unification to be realized. From 
the right panel in this figure, we see that the LHC has had a pretty significant effect on the 
expectation of M uni f. We clearly see that scales of order 10 8 — 10 15 GeV have become more 
favored, whereas there is a clear depletion of probability near the weak scale. 

Fig. shows the final product probability distribution for gf(M nni f) of the two distributions 
in Fig. [31 demanding that both expression in Eq. 13.91 agree. The result, not surprisingly, leads 
to a current preference towards small values of gf(M UVL if) or values of g\{M UTL n) > 0.4. 
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Figure 5: Final product probability distributions of the value of the gauge coupling g\ at the Gaugino 
Mass unification scale, demanding that both determinations of g\ shown in Fig. agree 
and apparent Gaugino Mass Unification takes place at the scale M un if. The number of 
bins is taken to be 50 for this plot, implying p(0\Flat) = 0.02. 
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4 General Gauge Mediation 



Gauge mediated SUSY-breaking is generically denned as a model in which supersymmetry 
breaking is transmitted to the observable sector via gauge interactions, leading therefore to 



flavor-blind parameters |60|-[65|. In Ref. [66| General Gauge Mediation (GGM) was defined 
as any theory in which all SUSY-breaking effects decouple from the MSSM in the limit of 
vanishing MSSM gauge couplings. 

The soft sfermion masses in GGM can be parameterized by a set of three parameters, Ai, 

3 

m} = Y,gtC i (f)A i , (4-17) 

i=l 

where the sum runs over the gauge groups of the MSSM. Here Ci(f) is the quadratic Casimir 
of the superfield / under the gauge group i, which, for a fundamental representation of SU(N) 
takes the value d(f) = (N 2 - l)/2N, while for U(l), C x {f) = Y 2 /A . Observe that we are 
implicitly working with a normalization of the gauge couplings consistent with their unification 
at the GUT scale, so g\ = 5gj(SM)/3, and Y 2 /4 = 3/5 x (Q - T 3 ) 2 . 

The gaugino masses are expressed in terms of three more parameters, Bi, given by 

Mi=g*Bi. (4.18) 

In order to generate a Higgsino mass parameter, //, and soft term, B^, of the correct order, 
gauge mediation may need to be supplemented by additional SUSY-breaking contributions in 
the Higgs sector. Therefore, we assume that in the case of the soft Higgs masses, the expression 
given in Eq. (I4.17P may be modified, 

m Hu = m ! 3 + 5 « ' 

m H d = m L 3 + 5 d. (4.19) 



Due to flavor independence, Db 13 and Dl 13 vanish in GGM. Moreover, the RGI D Xl also 
vanishes, as can be easily checked using its definition in Table HI The invariant Dz presents 
a simple dependence on the mass parameters and provides information on Sd- Therefore, the 
probability distribution for 5d in GGM can be read directly from the one of Dz presented before 
in Fig. HJ Of the other RGIs, there are six that probe the high-scale mass parameters of pure 
GGM, namely the J^s and the Im^- We shall mostly concentrate on those invariants in this 
section. 

As mentioned in Section 2.3, we exclude the other two invariants, Iy a and Dy 13H as they 
depend on too many parameters and currently it is difficult to obtain meaningful information 



1(3 



from them. Observe, that, eventually, the invariants Dy 1ZH and Iy a can be used to determine 
the gauge couplin gs a t the high scale and also probe possible non-universal corrections to the 



Higgs soft masses [39]. 



As emphasize above, we shall concentrate on the RGIs with explicit dependence on the gaug- 
ino mass parameters to extract information about the A{ and Bi. From the 1^ we immediately 
obtain 

Bi = I Bi , (4-20) 

and these distributions can be seen in Fig. [TJ 

In order to obtain information on the A t , both Ig. and Im, must be used: 

Al= 12f%_4 

33 V fl? 1 
M- 2(^-4 
2 i hi 3 T 2 



where the gi are the gauge couplings at the messenger scale, M. 

The probability distributions for the A\ for three different values of the messenger scale are 
given in Fig. EJ The shaded green area, as well as the green, red and black lines have the same 
interpretation as in Figs. [1] and [2j 

From Fig. [6j it is clear that positive and sizable values of A\ are preferred, independent 
of the messenger scale. On the other hand, values of A 2 close to zero are somewhat favored, 
although sizable values are equally likely. Finally, A 3 can be small and negative, but positive 
and sizable values are equally or more likely than the negative ones. One can check, that if 
the values of the A3 are the ones associated with the regions of maximal likelihood, which 
correspond to negative values of this parameter, and assume A 2 to be small, the boundary 
condition for the left-handed squarks square mass parameters would be negative, unless the 
largest values of A\ at each messenger scale are selected. These boundary conditions, together 
with the ones for the Bi, would then lead to somewhat light left-handed sfermions compared 
to the right-handed ones. 



5 Minimal Gauge Mediation 

Minimal Gauge Mediation (MGM) is a particular gauge mediated model in which the soft 
supersymmetry-breaking parameters are obtained through the interaction of messenger parti- 
cles that transform under the 5 + 5 representation of SU(5). The general assumption is that 
these messengers acquire mass via their interaction with a singlet superfield. When this su- 
perfield acquires a vev, it fixes the messenger scale, and its F-term, F$, fixes the scale of soft 



17 





p(f) | Flat Unweigh 

p(S | pre-LHC) 

h p(f) | poSI-LHC) 

----- p(8|Exp.) 

SF=0.21 


'h 




; i 


i 


\ 



x10° 



-30 -20 -10 



10 20 30 



p(B | Flat Unweighted)J~__L 




_ p(B | posl-LH 






p(B|Exp.) 

SF=0.23 






t t I i t t t I 1 1 


, , 1 , , , , 1 , , , , 1 , , , 


ll 



x10° 



A, (GeV) 2 withn= 10 b TeV 



-30 -20 -10 10 20 30 
A 3 (GeV) 2 with n = 1 8 TeV 



p{6 | Flat Unweight< 

- p(H | pre-LHC) 

- p[" | post-LHC) 

- p{9 I Exp.) 
SF=0.21 




x10° 



-30 -20 -10 10 20 30 
A, (GeV) 2 withn = 10 12 TeV 



Figure 6: Distribution of the GGM parameters A{ before and after the LHC constraints are added 
(green and red lines), flat distribution (shaded green) and subtracted probability distribu- 
tion (dashed black line). The three set of values are associated with different values of the 
messenger scale, M = 10 5 , 10 8 and M = 10 12 TeV. For the first and last row, no. bins = 50 
giving p(0\ Flat) = 0.02. 
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supersymmetry-breaking parameters through the identity 



M > = St ■ (5 ' 22) 

Hence in minimal gauge mediated models 

B < = B = iBs • (M3) 

with % — 1, 2, 3. Moreover, the scalar masses at the messenger scale are obtained at the two- loop 
level and acquire the value 



™;-£ 2 f^,(/)§f. (5.24) 

1=1 



Hence, in MGM 



2F 2 



Therefore, minimal gauge mediation is a model with just two parameters, B and the mes- 
senger scale, M (g 2 and g 3 can be written in terms of g± through the RGIs J S2 and I ga , which 
is then just a function of M). 

In order to determine the probability distribution for the parameters B and gf(M), we use 
the ones of the i^s and the J^s. From Table HJ we see that the I Mi must fulfill the following 
relations 



38gf(M)g 

J Mi 3 = U, 

5 

I M2 - 2g*(M)B 2 = 0, 

I M3 +2gl(M)B 2 = 0, (5.26) 

which together with the equations Jb 4 = i? define a system of 6 equations with only 2 unknowns. 
For every value of Ib { one can obtain a value of B that leads, from the 3 equations in Eq. f)5.26p . 
to 9 independent values of gf(M). In addition, these 3 equation lead to 3 different sets of 
simultaneous equations that can be solved for B and gf(M) independently. This leads to 
another 6 solutions for B and 3 for gf(M), leading to a total of 9 solutions for B and 12 for 

The set of 9 solutions that we use to compute the probability distributions for B are: 



19 



B x = I Bl , (5.27) 
B 2 = Ib 2 , (5.28) 
B 3 = I Bs , (5.29) 



IJi 5 — T F= 1 10. oil) 

V^(1089/ Ml - 95/m 2 ) 
6,7 ^(121/ Ml + 95J M3 ) ' 1 ' } 

The corresponding set of 12 solutions for the gauge coupling at the messenger scale, gf(M) 
are instead given by: 
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9u 
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5y— ~lMa 



(5.44) 



(llV2/ B2 - 5I g3 7=/^)' 



(5.45) 



12 



(llV2J Bs - 51,3^3 



The probabilities corresponding to each of these are plotted in Figs. IHlfTUl 
Observe that, as is apparent from Eq. (15.261) and Tabled MGM is associated with negative 
values of Im 3 and positive values of Im 13 ■ As can be seen from Fig. [TJ these values of the RGIs 
are not the most likely ones consistent with the present constraints. However, wether a given 
model is likely or not is a very scan dependent question and hence we will not address that 
here. Instead, the probability distributions for the MGM parameters are computed for those 
configurations for which these conditions are fulfilled. The final distribution is obtained by 
multiplication of the independent probabilities of the 9 Bi and 12 gf.(M) solutions given in 
Eqs. 15. 27115. 451 The results are depicted in Fig. [3 

The messenger scale may be obtained from the value of the gauge coupling at this scale 
by using Eq. (13.161) . replacing M uni f by M mess . However, in contrast to the gaugino mass 
unification scale, the messenger scale is always a physical scale and therefore expected to take 
values between tens of TeV and the GUT scale, or equivalently, gauge coupling values of 
0.2 < gj(M mess ) < 0.5. Fig. (JTj) shows that values of the gauge couplings <7i(M mess ) > 0.6 tend to 
be preferred, which lie outside the physical region. Considering only the physical range, values 
of the messenger scale close to the GUT scale are slightly preferred. The most probable values 
of the parameter B are about 1.25 TeV and 4.25 TeV. Using the relation Mj = gfB and the 
values of the gauge couplings at the weak scale, M mess ~ 1.25 TeV would lead to a bino mass 
of the order of 250 GeV, a wino mass of about 500 GeV and a gluino mass of about 1.5 TeV. 
The larger value of B would lead to gaugino masses 3.5 times heavier than these ones. 

6 Conclusions 

Supersymmetric extensions of the Standard Model provide a relationship between the weak 
scale and the scale of the supersymmetry-breaking parameters, rendering it stable under quan- 
tum corrections. In the MSSM, the SM-like Higgs particle is predicted to be light. The fast 
decoupling of the supersymmetric particles from the precision electroweak observables make 
the MSSM predictions consistent with those of the SM with a light Higgs, in full consistency 
with what current data seems to suggest. However, no direct hint of supersymmetric particles 
has been observed experimentally and hence no information of the structure and origin of the 
supersymmetry-breaking parameters is provided by current experiments, apart from perhaps 
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Figure 7: Final product distribution of the MGM parameters B and gf(M). 
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Figure 8: Distribution for the MGM parameter B before and after the LHC constraints are added 
(green and red lines), flat distribution (shaded green) and subtracted probability distri- 
bution (dashed black line). The different sets are associated with different probability 
distributions listed in Eqs. 15.271 
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Figure 9: Distribution of the gauge coupling, gf, at the messenger scale before and after the LHC 
constraints are added (green and red lines) , flat distribution (shaded green) and subtracted 
probability distribution (dashed black line) . The different sets are associated with different 
probability distributions given in Eqs. I5.45I 
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Figure 10: Distribution of the gauge coupling, g\, at the messenger scale before and after the LHC 
constraints are added (green and red lines), flat distribution (shaded green) and sub- 
tracted probability distribution (dashed black line). The different sets are associated 
with different probability distributions given in Eqs. 15.451 

the indirect hints provided by the anomalous magnetic moment and the Higgs mass range. 
Once additional information from direct searches becomes available, a method to determine 
the structure of supersymmetry-breaking parameters at the messenger scale, as well at the 
messenger scale itself would be desirable. RGIs provide such a method, establishing a direct 
relationship between the observables at the weak scale and the messenger scale parameters. 

In this article we have studied the probability distributions of a set of RGIs in the MSSM 
arising from symmetry arguments. The distributions are analyzed at the TeV scale by making 
use of the constraints coming from flavor physics, LEP and Tevatron searches, Higgs physics 
and the anomalous magnetic moment of the muon, and separately from those, by constraints 
provided by the LHC. We have used a flat prior for the soft supersymmetry-breaking masses, 
using a pMSSM approach. The current constraints already provide interesting features in the 
probability distributions. 

As an example of the application of the RGIs, we have used them to analyze the question 
of Gaugino Mass Unification and also the possible realization of General and Minimal Gauge 
Mediation. The methods described here are quite general and may be applied to analyze the 
ultraviolet properties of the MSSM parameters in other interesting supersymmetery-breaking 
scenarios. 

We noticed that the scale of Gaugino Mass Unification is not necessarily identified with the 
messenger scale, but it can provide non-trivial information on the realization of minimal models 
of supersymmetry breaking. GGM provides a well- motivated example of flavor independent, 
supersymmetry-breaking models. The probability distributions for the GGM parameters can 
be determined from those of the RGIs and present some interesting features as well. They also 
lead to information on possible non-universal Higgs mass parameters at the messenger scale. 
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The determination of the messenger scale in GGM through RGIs demands the measurement 
of both the first and third generation fermion masses as well as the Higgs masses, and hence 
it is not practical at this moment. We also analyze the more simplistic subset of models given 
by MGM. Since the entire model space of MGM is determined by only 2 parameters, we are 
able to extract information about the possible scale of SUSY particles as well as the messenger 
scale in this scenario. 

It is clear that although the analysis we describe already has interesting features, the prob- 
ability distributions of the RGIs will become particularly useful when the LHC starts revealing 
the presence of supersymmetric particles at the weak scale. In such a case, the probability 
distribution of the RGIs will become sharper and will start showing important features of the 
supersymmetry-breaking mass parameters at the messenger scale. Due to the higher cross sec- 
tions for the production of supersymmetric particles, the higher luminosities and the higher 
energy reach, the 8 TeV run this year will lead to relevant constraints on the supersymmetric 
particle masses. It could also lead to the first hint of the presence of supersymmetry, beyond 
the indirect ones associated with Higgs search results. It will be therefore very interesting to 
repeat the analysis of the RGI distributions once the 2012 results are available. 
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Appendices 

A Probability Re-weighting and Re-scaling 

We are interested in a quantity which quantitative reflects the probability distributions of gen- 
eral functions of the masses, given the probability distributions for the masses themselves. The 
Markov Chain Monte Carlo (MCMC) method is used to scan over the pMSSM parameters in 
the range considered to be probed at the LHC. For each point corresponding to a model, a 
likelihood is computed, given certain experimental constraints. Since the MCMC technique 
scans the given parameter space along the isocontours of likelihood due to preLHC constraints 
listed in Table [III the ratio of the number of points scanned for any given value of a parameter 
to the total number of points gives the probability for that parameter value. This probability 
for a given point can then be re-weighted by the postLHC likelihoods to compute the current 
probabilities. We note however that the boundaries defining the pMSSM region scanned, in- 
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troduce an artificial effect in the resulting probability distributions. In the following, we will 
describe a method that can be used for eliminating this effect. In this method, we make the 
assumption that the LHC (as well as the pre-LHC) measurements will not be able to shed any 
light on the pMSSM parameter regions that are not scanned due to kinematic constraints, and 
assign a flat probability to these insensitive regions outside the scan boundary. 

Let us consider a two dimensional probability distribution p(x,y\0) of parameters x and y 
defined in a box where the variables x and y vary in the ranges {x\, X2}, {y±, 2/2}; given some 
observables O. Assume that x and y have flat priors corresponding to the soft parameters that 
were scanned over in the MCMC We are then interested in the probability distribution of some 
function, 6(x,y), given O: p(8(x,y)\0). As explained in Section 2.2, the naive computation of 
this probability, especially using a flat prior pl(x,y) = constant, will heavily reflect the size 
of the box alongside any other inherent probability distribution of this function. The aim is 
to define a probability p(8(x,y)\0) such that, if there is no condition on 6(x,y), then a flat 
distribution is obtained for p (6(x, y)). Any variations of this flatness should be something that 
reflects the actual variation of the probability due to the effect of O rather than the effect of 
having a bounded box, as is the case in the example given in Section 12.11 

Let us assume that the box contains a bins in x and b bins in y. The flat distributions are 
defined such that in the absence of any additional condition: 



pl(x) = - (A.46) 
d(y) = I (A.47) 
vl&y) = ^ (A.48) 

pirn = E pfef) ( A - 49 ) 

All{x,y}:e(x,y)=6i 

Note that p{,(0) is defined as the distribution that would be obtained for 9(x, y) if x and y 
have flat priors. This is the distribution that is referred to in the text as "p(#|Flat Unweighted)". 
This distribution itself is generally not flat, but will have a distinct shape reflecting the boundary 
conditions of the original x, y variables. Analogously, the probability for 6(x, y) given O is 

p(0i\O)= £ p{x,y\0), (A.50) 
An{x, y y.e(x,y)=9i 

where this is referred to as "p(6*|pre/postLHC)" in the text. An easy way to normalize this 
probability to obtain a flat distribution for the function 8(x,y) in the absence of non-trivial 
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conditions is to weig ht each bin, S„ by l/pftfc): 

p * m0 ) oc . (A.51) 

The superscript A denotes the fact that this effectively gives the average probability per unique 
{x, y} combination for each However, this has the effect of washing out small effects on the 
probability distribution from O, when &i is such that a large number of unique combinations 
of {x, y} contribute to a given value of 8(x, y). 

We propose an alternative method. First, instead of taking the ratio we shall consider the 
difference: p(9i\0) —Po(0i)- Clearly this quantity is not always positive and cannot be identified 
with a probability distribution. It has, however, the property that it becomes positive whenever 
the probability of 9\ is enhanced by the observations O and negative in the opposite case. We 
shall, hence, define a renormalized distribution p R (9 i \0) in the following way 

p f (e m ) = Mas [p£(0O] (A.52) 
p R (9 l \0) oc p{9 i \O)+[p f {9 m )-p f {9 i )\, (A.53) 

which is always positive since p(9i\0) is positive and so is the quantity between brackets. The 
above quantity, Eq. (1A.53I) has a clear interpretation : Let us first stress that, by definition, 9 m 
is such that it has the largest number of unique combinations of {x, y} contributing to it, with 
x and y varying with a flat distribution in the box. Let's call k\ the number of combinations 
corresponding to 0j. Therefore, all 9i have a smaller number than 9 m , fc, < k m . This is 
reflecting the fact that for i ^ m, the range of the original variables scanned over, x and y, 
did not include all the combinations necessary to weight the % bin of 9 the same as m. We 
have made the argument that the values of x and y not scanned are ones that will not be 
affected by LHC measurements. Hence we propose that these combinations are given the same 
weight as p{{9i)/ki = 1/ab. This leads, after proper normalization, to nothing more than the 
last term, between square brackets, in Eq. flA.531) . and hence the quantity p R (9i\0) reflects the 
actual probability distribution of 9i given O, taking away the effect of the range of the original 
scan. For this quantity to represent a probability distribution in the strict sense, it must be 
normalized to 1. Since p(9i\0) and p Q {9i) are quantities which are normalized to 1, assuming 
that the function 9i is evaluated in / different bins, the normalization factor is nothing more 
than C = l/(pl(9 m )l). Hence the properly normalized probability distribution for 9 is given 
by: 

^1°) = -nTTTl 1^1°) + \dVm) - dVi)] } ■ (A.54) 

P {9m) I 1 1 ]) 

We can see that this behaves the way we expect it to, by noting that when O has not impacted 
the probability of 9, i.e. p(9i\0) = pl(9i), p R (9i\0) = 1/7, so we obtain a flat distribution. On 
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the other hand, if the Po(#i) * s a constant, meaning that 0, has a flat distribution in the same 
flat basis as the original variables x and y, then pl(9i) = pl(9 m ) = 1/1 and we recovers p(9i\0) 
without any modification, as we should. 

In order to emphasize the impact of the experimental constraints in a more clear way, 
however, we have gone a step further. Since we assumed that the probability outside the range 
we scanned is flat, the ratio of the difference of any two probabilities from flat, (p R (9i\0) — 
l/l)/(p R (9j\0) — 1//), will remain invariant if we extended the range of the original scan, 
increasing the box size. Therefore, this quantity is than also invariant under an overall rescaling 
of the differences with the flat probability. 

Let us assume that there is a non-trivial impact of experiments on the RGI distributions, 
namely p R (9i\0) ^ l/l for at least one %. Considering 

p R (9 n \0) = Min [p R (9i\0)] (A.55) 

we define a scale factor, SF, such that the difference of this minimum with l/l is scaled to be 
l/l: 

SF^-p R (9 n \0^ = j (A.56) 
SF" 1 = l-p R (9 n \0)l. (A.57) 

We use the scale factor above to define a modified distribution 

1" 



p ss (9i\0) = \ + SF 



(A.58) 



I 

= -HTTT, \po^ + SF \pM°) ~ dVi)] } ■ (A-59) 
P {9 m )l L L JJ 

Once the scale factor SF is given, it is easy to translate this modified distribution, Eq. ( 1A.59j) to 
the original one, Eq. (1A.54I) . The quantity p ss (9i\0) has the virtue that when for a particular 
bin p R (9i\0) = l/l, meaning O has had no impact on the 9i probability, one obtains p ss (9i\0) = 
l/l. On the other hand when p R {9i\0) = p R (9 n \0), meaning when O has maximally decreased 
the probability for that 9 U p ss (9i\0) = 0. 

The fact that p ss (9i\0) will be invariant under a change in scan range of the original 
variables can be seen by inspecting Eq. IA.58I and noting that under a change of scan range, 
p ss (9i\0) = l/l when p R (9i\0) = l/l and by definition p ss (9 n \0) = 0. 

Even though p ss {9i\0) cannot be technically defined as a probability, it quantitatively re- 
flects the actual impact of O on the probability distribution of 9 in a way which is independent 
of the artificial impact of scanning a finite region, and, as stressed above may be easily con- 
nected with p R (9i\0), Eq. f[P4l) . We ran extensive numerical checks to make sure that this 
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quantity indeed behaves in the expected manner. We have therefore used p ss (6i\0) to repre- 
sent the probability distribution of the RGIs, giving the associated scale factor SF for every 
RGI distribution. In the text, in order to be more explicit about the meaning of these dis- 
tributions, p{9\0) was renamed "p(6*|pre/post-LHC)", while p (0\O) was renamed u p{9\Exp) 
Reweighted" . 



B pMSSM Parametrization 



The pMSSM, a 19-dimensional realization 67| of the R-parity conserving MSSM with param- 



eters defined at the SUSY scale, M SU sy = ^ m t 1 m t 2 i employs only a few plausible assumptions 
motivated by experiment: there are no new CP phases, the sfermion mass matrices and trilinear 
couplings are flavor-diagonal, the first two generations of sfermions are degenerate and their 
trilinear couplings are negligible. In addition, we assume that the lightest supersymmetric par- 
ticle (LSP) is the lightest neutralino, x?- We thus arrive at a proxy for the MSSM characterized 
by 19 real, weak-scale, SUSY Lagrangian parameters: 

• 3 gaugino mass parameters Mi, M 2 , and M 3 ; 

• the ratio of the Higgs vevs, tan = v^/vx] 

• the higgsino mass parameter, /i, and the pseudo-scalar Higgs mass, m^; 



• 10 sfermion mass parameters mp, where F = Q\, Ui, Di, L\, E\, Q3, U3, D 3 , L 3 , E 3 
(imposing m ( g l = m^, = mi 2 , etc.); and 

• 3 trilinear couplings A t , Ab and A T , 
in addition to the SM parameters. 



For each pMSSM point, Sof tSUSY3 . 1 . 6 [68] was used to compute the SUSY spectrum, 



Superlsov3.0 [69J was used to compute the low-energy constraints, micrOMEGAs2.4 [51] was 
used for the SUSY mass limits, and HiggsBounds2 . . 50] for the limit on the h° mass 3- 



Moreover, SUSYHIT (SDECAY1.3b, HDECAY3 . 4) [70] was used to produce SUSY and Higgs 



decay tables, and micrOMEGAs2.4 51] to compute the LSP relic density and direct detection 



cross sections. The various codes were interfaced using the SUSY Les Houches Accord 



C Soft Mass Parameters and RGIs 

As mentioned in Section 2, one can make use of the RGIs and three independent masses to 
determine all other soft breaking masses. As an example, we write down 2 sets of solutions 

1 In evaluating the Higgs mass limit, a Gauss-distributed theoretical uncertainty of a = 1.5 GeV was applied 
to the rrih computed with SoftSUSY, cf. row 8 in Table ITT1 and row 12 in Table Hill 
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with different unknown masses. All the masses and gauge couplings are at the same scale. The 
gaugino masses in both cases are given by 



Mi = I B .g* 2 = 1,2,3. (C.60) 
We write the first set of solutions in terms of 3 third generation masses: mg 3 , m u , 3 and m e . 3 , 



m 2 _ D B ri _ Dz_ _ 5I A/l 37m 2 4Jm 3 _ _ 247£>y 13g 

^ 2 2 2 66 2 3 Ll3 220 

D X1 3I Ya gl 5 2 4 3 2 4 4 2 4 3m 2 3 

+ 40 + 22 + 66 2 2 3^3^ + 2 , 



(C.61) 



2 _ 3ggig -Dz 27mi 4Jm 3 £> Ll3 13£>y 13H 3£> Xl 

m ^ " 2 + ^3~ " ^ 2 + ~3 2 44~~ + ~T 

" ^ " ^ + ^ " 3^ + ^ + 3< - ^ , 



(C.62) 



«3 



^ 13 + ^-3/ M2 -^^ + ^-^!--/l^ + 34^ + 2m 2 -m 2 , 

(C.63) 



m Qi = 3^ (20J Ml + 5940/ A/2 - 3520J M3 + 78D Yl3H - 627D Xl + 60Jy a ^ - 20/^^ (ag4) 



- 5940/i 2 £ 2 4 + 3520/I 3 ^ 4 ) 



™£ 3 = ^ {-Wh h + 330/m 2 - 110D Ll3 - 26D Yl3H - 11D X1 - 20I Ya g 2 + lOI 2 ^ 4 (c 65) 



m 



1 



440 



- 330/^ 2 4 + 110my 



(20J Ml + 660/ A/2 - 26D Yl3H - UD X1 - 20gl (l Ya + I 2 Bi gl) - 660/l^ 4 ) , (C.66) 



m 



dl 1980 



(40/ Ml - 1760/ M3 + 78£>y 13H + 33L> X1 + 60I Ya g* - 40I 2 Bi gt + 1760/J^ 4 ) , (C.67) 



ml = (80I Ml - 880/m 3 - 78D Yl3H - 33D Xl - 60/y^ 2 - 801^ + 880/| 3 ^) , (C.68) 



990 
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Alternatively, the second set of solutions is given in terms of the 2 soft masses for the Higgs, 
m Bu and m^, and a third generation squark mass, rriQ 3 : 



2 D Bl3 D z 5J Ml 8/ M3 2D Ll3 247Dy 13H 

* * 99 M2 9 3 330 (C.70) 



60 11 99^1+^2 + 9^3+ 3 : 



ie 2 = -4D Bl3 + 2D Z + ^ + 3J M2 _ 1^3 + 3 /J Ll3 + _ 1^ 

33 2 3 13 55 5 (C.71) 



2 V#1 1 r2 4 o r 2 „4 , 16 r2 4 



+ Z1 Tf ± ~ 55 r ^i - 3/ ^2 + + 2m| d + 2m^ - 6m 2 



11 33 

1 



Qi = ^ (20V + 5940/ M2 - 3520V + 78D Yl3H - 627D Xl + 60I Ya gl - 20/^ (c ?2) 

- 5940I 2 2 <7 2 4 + 3520/ 2 3 <? 3 4 ) , 



2 _4D Bl3 D z 4I Ml 
»»*--3 3" + ^9"" 2/M2 + 

+ I -^--I 2 B gt + 2I 2 B g i 2 - 
33 99 



m L, — ~ ZJJ B 13 + i->Z ^77" + <J-*M 2 o r i/L 13 



8V 3 


2D Ll3 


91Dy 13H 


19D Xl 
1 — 


9 


3 


110 


60 


-/ 2 o 4 

g B^iJZ 


2m 2 „ 
i« 

3 


^2m 2 Q3 , 





(C.73) 



m 



2 

Ll ~ 440 



20 33 

1 



33 2 3 " 10 (C.74) 

QD 8 

yiy Xi i- L r2 4 or2 4, r2 4, 2 , 2 o2 

V^i - 3/ B2 £ 2 + -/ B3 ^ 3 + m H d + m H u ~ 3m Q3 , 



(20V + 660V - 26D Yl3H - 11D X1 - 2Qg\ (l Ya + I 2 Bl gf) - Q60I 2 B2 g%) , (C.75) 



m 



* 1 
rfl 1980 



(40V " 1760V + 78Dy 13ff + 33D Xl + 60V<? 2 - 40/J^ 4 + 1760/ 2 3 ^ 4 ) , (C.76) 



m 2 



Ul g^ (80V " 880V " 78D Fl3H - 33D Xl - 6QI Ya £ ~ 80/f^ 4 + 880/ 2 3 s 4 ) , (C.77) 
< = ( 40/Ml + 26Dy -« + UD ^ + 20Jy ^ ~ 40J ^^) (C - 78) 
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D Flavor-Blind Models 



The most immediate consequence of flavor-blindness is the vanishing of D # 13 and -D_l 13 . There- 
fore these invariants provide us with a direct test of the flavor-independent hypothesis with 
a minimal set of measurements. More precisely, they allow this hypothesis to be ruled out: 
measuring Db 13 7^ or D r„ ^ at the low scale implies high-scale family non-universality; 
however, as noted in Ref. [33j, measuring Db 13 = and Z?z 13 =0 at the low scale does not 
necessarily indicate high-scale universality. 

Current experimental data from flavor physics strongly motivates a flavor- universal medi- 
ation mechanism for SUSY-breaking. Accordingly, if Db 13 and Dl 13 are found to vanish, it is 
reasonable to proceed a step further and attempt to extract constraints on the high-scale values 
of the flavor-blind MSSM soft parameters from the RGIs. 

The 7 scalar and 3 gaugino soft mass parameters in the flavor-blind MSSM can be expressed 
uniquely in terms of the 10 invariants D Xl through Im z listed in Table [B These are listed in 
Eqs. flC.60p and flD.79p - flD.85p . Note that these relations depend on the 3 gauge couplings and 
further all couplings and soft parameters are assumed to be given at the messenger scale: 



m\ = 



^(^D Yl3H + 11D X1 + 20(^4 + 33g 4 2 I B2 ) - (h h + 33/ M2 ) + y?Iy«)) , (D.79) 

m k =<- \ D z , (D.80) 

m k = m L- \ D z - n^3H + , (D.81) 

ml = ^( 26D ri*x + 11D x. ~ 20(2(^4 - I Ml ) - g 2 I Ya )) , (D.82) 

^(78LV 13H + 33D X1 + 20(4(^4, - HflgiSJ - (h h - H/a/J) + 3g 2 I Ya )) , 

(D.83) 



2 1 

m u = —■ 



m\ = Y^(78D YlSH + 33D X1 - 20(2(^4 - Ug 4 3 I 2 B3 ) ~ (h h - 441*,)) - ?>gll Ya )) , 

(D.84) 

m% = ——(78^^-627^. 
Q 1 3960 v 13H Xl 

-20((^ + 297g 4 2 I 2 B2 - 176^ 4 /1 3 ) - (I Ml + 297/ A / 2 - 176/^) - Sg 2 I Ya )) .(D.85) 

Using the invariants I g2 and I g3 these may be expressed entirely in terms of g±. Equivalently, 
one can reduce the degrees of freedom at the high scale to a single parameter, which can be taken 
to be the value of that scale. In particular this permits tests of more restrictive flavor-universal 
models such as mSUGRA, taking g\ at the GUT scale. 
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